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ABSTRACT 

O This paper reports the latest results from a millimeter-wave (CO) spectroscopic survey 

of IRAS-detected radio galaxies with Li.4GHz ~ i o 23 ~ 28 W Hz -1 in the redshift range 
z ~ 0.02 — 0.15. The IRAS flux-limited sample contains 33 radio galaxies with different 
radio morphologies and a broad range of infrared luminosities (Ltr = 10 9-12 L ), allowing 
for an investigation of (a) whether low- z radio-selected AGN reside in molecular gas-rich 
^ . host galaxes, and (b) whether the CO properties are correlated with the properties of 

CN | the host galaxy or the AGN. All of the radio galaxies in Mazzarella et al. (1993) and 

Mirabel et al. (1989) have been reobserved. Three new CO detections have been made, 
raising the total number of CO detections to nine and setting the survey detection rate 
at ~ 25%. Many of the CO lines have double-peaked profiles, and the CO line widths are 
broad (average At> FWH M ~ 500 ± 130 km s _1 ), exceeding the average CO widths of both 
ultraluminous infrared galaxies (300±90 km s _1 ) and Palomar-Green QSOs (260±160 km 
s" 1 ), and thus being indicative of massive host galaxies. The CO luminosities translate 
into molecular gas masses of ~ 0.4 — 7 x 10 9 M , however, the 3cr CO upper limits for 
nondetections do not rule out a molecular gas mass as high as that of the Milky Way 
(~ 3 x 10 9 M ). Optical images of eight out of nine molecular gas-rich radio galaxies 
show evidence of close companions and/or tidal features. Finally, there is no obvious 
correlation between radio power and molecular gas mass. However, it is notable that 
only one F-R II galaxy out of 12 is detected in this CO survey; the remaining detections 
are of galaxies hosting F-R I and compact radio jets. 

. ^ [ Subject headings: galaxies: active — galaxies: interacting — galaxies: ISM — ISM: 

^ ■ molecules — infrared: galaxies - 
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wavelength surveys of luminous infrared-selected 
galaxies, which have shown these galaxies to be 
merging/interacting disk galaxies and have shown 
an increasing fraction of infrared galaxies with 
optical, near-infrared, and X-ray signatures of 
imbedded AGN as a function of increasing lu- 
minosity (e.g., Sanders et al. 1988a,b; Veilleux, 
Sanders, & Kim 1997, 1999; Veilleux et al 1999; 
Vignati et al. 1999). 

Star formation appears to be a component of 
the activity occuring in the host galaxies of AGN 
- a growing number of radio galaxies and QSOs 
have detectable amounts of molecular gas and 
young UV-luminous knots indicative of recent 
star formation (Philips et al. 1987; Sanders, 
Scoville, & Soifer 1988c; Lazareff et al. 1989; 
Mirabel, Sanders, & Kazes 1989; Mazzarella et 
al. 1993 - hereafter, Paper 1; Reuter et al. 1993; 
Evans et al. 2001; Allen et al. 2002; Scoville 
et al. 2003). In addition, the fraction of classi- 
cal AGN with dust in their host galaxies is sub- 
stantial (e.g., Golombek et al. 1988; Sanders et 
al 1989), with ~ 50% of 3C radio galaxies and 
~ 70% of Palomar-Green (PG: Schmidt & Green 
1983) QSOs at z ^ 0.12 having IRAS 8 detections 
at 60/iin. The Infrared Space Observatory (ISO) 
detected thermal dust emission in a larger frac- 
tion of radio galaxies and QSOs (Meisenheimer 
et al. 2001; Haas et al. 2003, 2004), providing 
further proof of rich ISMs in the host galaxies of 
radio galaxies and QSOs. This dust is most likely 
a combination of dust intrinsic to the progenitor 
galaxies, as well as dust produced through stellar 
mass loss from a recent /ongoing starburst. 

The importance of star formation in AGN host 
galaxies is also inferred via recent black hole sur- 
veys in nearby galaxies. These observations have 
shown that the massive spheroidal component of 
many nearby normal galaxies contain quiescent 
supermassive nuclear black holes (M. = 10 6 ~ 9 
M : e.g. Magorrian et al. 1998). A key re- 
sult that has come from these surveys is the scal- 
ing of black hole mass with stellar bulge mass 
(e.g., Magorrian et al. 1998). Such a correla- 
tion may indicate that the formation of stars in 
the bulge component of galaxies is intimately as- 
sociated with the building of supermassive black 
holes. 

While the molecular gas properties of infrared- 
selected galaxies have been studied in detail, CO 

8 I.e., the Infrared Astronomical Satellite. 



observations of large samples of radio galaxies 
and QSOs have lagged behind. Thus, conclu- 
sions about the CO properties of classical AGN 
and their host galaxies is often based on small 
number statistics or incomplete data sets. The 
present study focusses on the global molecular 
gas properties of a sample of low-redshift radio 
galaxies, and it serves as an extension of the sur- 
vey presented in Paper 1. The primary issues 
to be addressed with this survey are (a) whether 
low- ,2 radio-selected AGN reside in molecular gas- 
rich host galaxies and (b) whether the CO prop- 
erties are correlated with the properties of the 
host galaxy or the AGN. To address the latter 
issue, the millimeter data are complemented by 
infrared, radio, and X-ray data to allow a com- 
parison between the available amount of fuel for 
star formation (and perhaps AGN activity) and 
the energy emitted from the AGN and/or host 
galaxy at these wavelengths. The analysis pre- 
sented here builds primarily upon similar analy- 
sis done of low-redshift , infrared-selected galaxies 
over the last twenty years (e.g., Sanders & Mirabel 
1996 and references therein; Franceschini et al. 
2003; Trentham et al. 2005) and recent analysis 
of a sample of CO-detected PG QSO host galax- 
ies (Evans et al. 2001). These galaxies, which 
may emanate from the same parent population 
as the radio galaxies or be linked via evolution, 
will serve as the primary comparison samples. In 
addition, optical imaging of a radio galaxy sub- 
sample will be briefly discussed. Parts of this 
extended millimeter survey have been presented 
elsewhere (Evans 1996; Sanders & Mirabel 1996; 
Evans 1998; Evans et al. 1999a,b). 

This paper is divided into six sections. The se- 
lection criteria of the radio galaxy sample is dis- 
cussed in §2. Section 3 is a summary of the obser- 
vations and data reduction. In §4, the CO(l — > 0) 
emission-line properties are presented, along with 
the method of calculating the molecular gas mass. 
Section 5 contains a brief summary of the opti- 
cal imaging data and an extended discussion of 
the infrared, radio, CO, and X-ray properties of 
the radio galaxies relative to those of infrared- 
selected galaxies, PG QSOs, and elliptical galax- 
ies. Section 6 is a summary. Throughout the 
paper, H = 75 km s _1 Mpc -1 , q = 0.5, and 
A = 0.0 are assumed. 

2. SAMPLE 



EVANS ET AL. 



3 



The radio galaxy sample was selected primar- 
ily from a list of radio galaxies correlated with 
the IRAS archive by Golombek, Miley & Neuge- 
bauer (1988). The radio galaxies were selected 
for CO observations based on the following crite- 
ria: (i) a declination, 5 > —20°, to allow observa- 
tions from both the Kitt Peak 12m Telescope and 
the James Clerk Maxwell Telescope (JCMT), (ii) 
IRAS flux densities at 60/im, fm^m, or 100/xm, 
Zioo^m, greater than 0.3 Jy, with the assumption 
that an IRAS detection at this level is a good indi- 
cator of a dusty, gas-rich interstellar medium, and 
(Hi) a redshift in the range 0.02 ^ z <; 0.15, to en- 
sure detectability of the molecular gas mass of the 
Milky Way (~ 3 x 10 9 M ), if present, with single- 
dish millimeter telescopes used for the survey 
(for the sensitivity calculation, the assumption is 
made that the CO lines have a full width at half 
the maximum intensity width of Awfwhm ~ 250 
km s _1 , which is the average value for a sample 
of IRAS galaxies observed by Sanders, Scoville, & 
Soifer 1991). Such criteria yield 31 radio galaxies. 
Two additional radio galaxies, PKS 0624-20 and 
3C 285, were identified via a search of the NASA 
Extragalactic Database (NED) to have properties 
consistent with the selection criteria; their ad- 
dition to the sample raises the total number of 
galaxies to 33. Note that two of the 33 galaxies, 
the relatively nearby galaxies NGC 6251 and 3C 
272.1, have CO non-detections in published sur- 
veys by Elfhag et al. (1996) and Knapp & Rupen 
(1996), and thus were not reobserved as part of 
this survey. General properties of the sample of 
33 radio galaxies are provided in Table 1. 

3. OBSERVATIONS AND DATA REDUCTION 



3.1. CO Spectroscopy 

Millimeter observations of the redshifted 
CO(l — > 0) emission in 31 of the 33 radio galaxies 
were made with the NRAO 9 12m Telescope dur- 
ing six observing periods between 1992 Decem- 
ber and 2000 June. The telescope was configured 
with two 256 x 2 MHz channel filterbanks and 
dual polarization SIS spectral-line receivers tuned 
to the frequency corresponding to the redshifted 
CO(l — > 0) emission line. (Optical emission line 
redshifts were adopted.) Observations were ob- 
tained using a nutating subreflector with a chop 
rate of ~ 1.25 Hz. Six minute scans were taken, 
and a calibration was done every other scan. Dur- 

9 The NRAO is a facility of the National Science Foundat 
Universities, Inc. 



ing the course of each observation, data were also 
taken with the bandpass centered at velocity off- 
sets of ± 20 or 50 km s _1 of the velocity corre- 
sponding to the redshift of the optical emission 
lines. Shifting the velocities in such a manner 
minimizes any ripples inherent in the baseline. 
Pointing was checked every few hours by observ- 
ing standard continuum sources and was generally 
determined to be accurate to within a few arcsec- 
onds. 

Observations of the redshifted CO(l — > 0) emis- 
sion in Cygnus A and PKS 1345+12 were done 
at the IRAM 30m telescope during 1995 Decem- 
ber and 2003 December observing periods, re- 
spectively. The 3 and 1 millimeter receivers were 
used simultaneously in combination with the fil- 
terbanks and an autocorrelator to provide band- 
widths of 512 and 600 MHz, respectively. Dur- 
ing observations, the pointing was monitored by 
observing the planets and standard continuum 
sources. 

Observations of the redshifted HCN(1 — > 0) 
emission in PKS 0502-10 and B2 0722+30 were 
done at the IRAM 30m telescope during the 2003 
December observing period. The setup was iden- 
tical to the CO(l — > 0) observation of PKS 
1345+12 obtained during the same period. HCN 
observations of the strong 3mm continuum galax- 
ies 3C 84 and 3C 120 were also attempted, but 
have been left out of the present discussion due to 
the presence of large baseline ripples in the data. 

Observations of the redshifted CO (2 — > 1) emis- 
sion in PKS 0502-10 and B2 0722+30 were done 
with the James Clerk Maxwell Telescope (JCMT) 
in 1995 September. The spectral line receiver 
(A2) were used together with the Digital Auto- 
correlation Spectrometer in wide-band mode (750 
MHz bandwidth). Because of excess noise near 
the edge of the passband, the usable bandwidth 
was only ~ 700 MHz, which corresponds to a total 
velocity coverage of ~ 900 km s _1 at 230 GHz. All 
observations were obtained using a nutating sub- 
reflector with a chop rate of ~ 1.25 Hz. Data were 
stored as six-minute scans, and a chopped wheel 
calibration was performed after every other scan. 
Pointing was monitored every few hours, and ve- 
locity shifts of ± 20 km s" 1 were made to min- 
imize baseline ripples. A journal of observations 
for all of the single-dish observations are provided 

ion operated under cooperative agreement by Associated 
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in Table 2. 

The data reduction for all of the single-dish mil- 
limeter and submillimeter observations were re- 
duced using the IRAM data reduction package 
CLASS. The individual scans were checked for 
baseline instabilities, then averaged together. A 
linear baseline was subtracted from each spec- 
trum; the cases where emission lines were de- 
tected, the baseline subtraction was done at ve- 
locities outside of the range of the emission line. 
Each spectrum was smoothed to ~ 20—50 km s _1 . 
Finally, line fluxes were measured by numerically 
integrating over the channels in the line profile, 
and the line widths were measureed as full width 
at 50% of the peak flux. The resultant spectra 
are plotted in Figures 1-3. 

3.2. Optical Imaging 

Ground-based imaging observations of 7 of the 
radio galaxies in the sample were made at the UH 
2.2m Telescope. The U', B, and I-band images 
were obtained over the period 1998 March 25 to 
2000 December 31 using the Orbit Semiconduc- 
tor 2048x2048 CCD camera. The original scale 
is 0.09'Vpixels (Wainscoat 1996), but the CCD 
was read out with 2x2 pixel binning. Three to 
four dithered exposures were taken, with integra- 
tion times ranging from 360 to 720 seconds each 
for U', B, and I-band filters. 

The U', B, and I-band data reduction was per- 
formed using IRAF. The data reduction consisted 
of flatfielding individual images, scaling each im- 
age to its median value to correct for offsets in 
individual images, then shifting and median com- 
bining the images. The final images, boxcar 
smoothed by 4 x 4 pixels, are shown in Figure 
4. 

4. RESULTS 

Table 3 is a summary of the CO emission-line 
properties of the 31 radio galaxies observed as 
part of this survey. All of the radio galaxies ob- 
served for Paper 1 were reobserved, as well as two 
additional radio galaxies previously detected in 
CO (3C 84: Lazareff et al. 1989; Mirabel et al. 
1989; Reuter et al. 1993, PKS 1345+12: Mirabel 
et al. 1989). In one case (TXS 1506+345), the 
pointing was shifted by ~ 30" from that used in 
Paper 1 in order to center the galaxy in the beam. 
Nine radio galaxies are detected in CO(l — > 0), 

10 Note that a detection of 3C 31 has also been reported 



including three new detections (3C 31, PKS 0502- 
10, 3C 293) 10 . However, the previously reported 
detection of CO emission in B2 0648+27 was not 
confirmed. The radio galaxy PKS 1345+12 was 
the weakest CO detection in the sample, thus in- 
dependent CO observations of the galaxy were 
also obtained with the IRAM 30m Telescope to 
improve the determination of the CO flux and line 
shape (Figure lb). The average full width at half 
the maximum intensity velocity, AfpwHM, of the 
detected galaxies is 500 km s" 1 . Several of the 
galaxies have double-peaked profiles. 

Table 3 lists the measured CO intensities and 
luminosities. The CO luminosities of 3C 120, 
B2 0722+30, B2 1318+34, and TXS 1506+345 
agree with the values in Paper 1 to within ~ 
20%, and the values of Aw FWH m for all except 
TXS 1506+345 are similar (TXS 1506+345 is re- 
ported as having a 2-component CO emission line 
profile in Paper 1). There is a major discrepancy 
between the new CO measurements of Perseus A 
and PKS 1345+12 and those in the literature. For 
Perseus, the line width of CO luminosity are in 
agreement with the measurements of Lazareff et 
al. 1989 and Reuter et al. (1993), but L' co is 
a factor of two higher than the value determined 
by Mirabel et al. (1989). The discrepancy with 
the Mirabel et al. (1989) measurements results 
from what appears to be extended emission in the 
wings of the line in Figure 1 - if this emission is ig- 
nored, then the CO luminosity is 6 x 10 8 K km s _1 
pc 2 , within 22% of the value reported by Mirabel 
et al. (1989). The PKS 1345+12 CO luminosity 
comparison is more difficult to rectify - the CO 
luminosity is a factor of two lower than measured 
by Mirabel et al. (1989). Note, however, that the 
errors associated with the Mirabel et al. (1989) 
data are not known. The present PKS 1345+12 
CO measurements are correct to within the er- 
ror of those reported in Evans et al. (1999) if a 
20% calibrated error on the present observations 
is assumed. 

Many of the observations resulted in non- 
detections. For each of these galaxies, upper lim- 
its to the CO intensity were calculated via, 



3T rms Aw FWH M Uz n -li /-.n, 

T mh Av < ; [K km s J, (1; 



^AfFWHM/A? 



where T mb is the main beam temperature, 
by Lim et al. (2000). 
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AfpwHM is the average FWHM velocity for the 
detected radio galaxies (= 500 km s -1 ), and T rms 
is the root-mean-squared main beam temperature 
of the spectral data for a velocity resolution of 
Av Tes . 

To calculate the CO luminosity, the luminos- 
ity distance for a source at a given redshift and 
a A = Universe was first calculated using the 
equation 



D L = cH^q^ 2 \zq Q + (q - 1) (^2q z + 1 - l) } 

[Mpc]. (2) 

Given the measured CO flux, Sco^v [Jy km s _1 ], 
the CO luminosity of a source at redshift z is 



CO 



S co AvDl(l + zy 



[K km s 1 pc 2 ], 



(3) 



(Solomon, Downes, & Radford 1992) where c [km 
s _1 ] is the speed of light, k [J K" 1 ] is the Boltz- 
mann constant, and z/ obs [Hz] is the observed fre- 
quency. In terms of useful units, -^co(i->-o) can ^ e 
written as 



L- co = 2.4xl0'f SC ° A \)(£ 

co v j y km s / \ m p c J 



(l + z) 



-1 



[K km s" 1 pc 2 ]. 



(4) 



To estimate the mass of molecular gas in these 
radio galaxies, the reasonable assumption that the 
CO emission is optically thick and thermalized, 
and that it originates in gravitationally bound 
molecular clouds, is made. Thus, the ratio of the 
H 2 mass and the CO luminosity is given by 



a 



M(H 2 



(X 



n(H 2 ) 



[M (K km s^ 1 pc 2 



(5) 



where n(H 2 ) and T b are the density of H 2 and 
brightness temperature for the CO(l— >0) transi- 
tion (Scoville & Sanders 1987; Solomon, Downes, 
& Radford 1992). Multitransition CO surveys of 
molecular clouds in the Milky Way (e.g. Sanders 
et al. 1993), and in nearby starburst galaxies 



(e.g. Glisten et al. 1993) have shown that hot- 
ter clouds tend to be denser such that the den- 
sity and temperature dependencies cancel each 
other. The variation in the value of a is approx- 
imately a factor of 2 for a wide range of kinetic 
temperatures, gas densities, and CO abundance 
(e.g. a = 2 - 5M Q [K km s" 1 pc 2 ]" 1 : Rad- 
ford, Solomon, & Downes 1991). Solomon et al. 
(1997) and Downes & Solomon (1998) have made 
use of dynamical mass estimates of a low-redshift 
infrared galaxy sample observed in CO with the 
Plateau de Bure Interferometer to argue that a 
may, in some cases, be as low as 1 M Q (K km s" 1 
pc 2 )" 1 . 

For the present sample of galaxies, an a = 
1.5M (K km s _1 pc 2 ) -1 is adopted, which is the 
approximate value derived for luminous and ultra- 
luminous infrared galaxies (Solomon et al. 1997; 
Downes & Solomon 1998; Evans et al. 2002). 
However, note that a for some, or all, of these 
galaxies may be a factor of 2 — 3 higher (i.e., the 
approximate value derived for bulk of the molec- 
ular gas in the disk of the Milky Way: Scoville & 
Sanders 1987; Strong et al. 1988). The derived 
molecular gas masses are listed in Table 3. Taking 
these masses and comparing them with the dust 
masses, Md us t, which are calculated via 



M, 



dust 



4.78 



/ /lOO/im \ / Dl 



= 143.88/T duat 



V Jy J \ y Mpc / 



[Ms] (6) 

(see Young et al. 1989) where the dust tempera- 
ture, Td U st is calculated from the 60/im and 100/im 
flux densities via 



T, 



dust 



-a+z) 



82 



ln(0.3/60^m/ /l00//m) 



-0.5 



[K], 



(7) 



yields an average molecular gas-to-dust mass ra- 
tio of 860 ± 460 (see Table 4). By comparison, 
the gas-to-dust ratio is ~ 600 (a = 4.8M [K km 
s _1 pc 2 ] -1 ) for IRAS-detected spiral and luminous 
infrared galaxies (Young & Scoville 1991 and ref- 
erences therein). Note that the estimated dust 
masses in Table 4 are likely lower limits; the bulk 
of the dust in these galaxies may be in a colder 
(< 30K) dust component (Young & Scoville 1991; 
Devereux & Young 1990). 



6 



CO IN RADIO GALAXIES 



Given the adopted value of a, it is worthwhile 
to make a conservative estimate of the dynami- 
cal mass in order to determine the percentage of 
the total mass in the form of molecular gas. The 
dynamical mass within the region of the galaxies 
containing molecular gas and dust can be calcu- 
lated via 



Mdyn ~ A "™ fico = 226 ^ Au ™ M V 



G 



km 



-i 




[M ], 



where A-u F whm is the full CO velocity width at 
half the maximum flux density (Table 3) and Rco 
is the radius of the CO distribution. The quan- 
tity Rco is estimated by assuming optically thick, 
thermalized gas with a unity filling factor and a 
blackbody temperature, Tbb, equal to that of the 
dust. Thus, the radius of the CO distribution is 
calculated via 



-Rco 



u 



CO 



TrTbbA^FWHM 



[pc] 



(9) 



(see Table 4). The resultant size estimates 
yield M dyn ~ 0.22 - 2.5 x 10 10 M and thus 
M(H 2 )/M dyn = 0.10 - 0.94 with a mean value 
of 0.45 ±0.29. 

Note that the values of Rqo are likely to be 
underestimated, and thus the derived values of 
M(H 2 )/M dyn are lower limits. For example, the 
column densities derived from Rco and M(H 2 ) 
are in the range 1 — 4 x 10 24 cm -2 . In contrast, 
the X-ray column densities estimated from X-ray 
observations of 3C 120 and PKS 1345+12 are in 
the range of 0.2 — 4 x 10 22 cm" 2 (Sambruna et 
al. 1999; O'dea et al. 2000). Further, the CO 
disks in 3C 293 (Evans et al. 1999a) and 3C 31 
(Okuda et al. 2005) are extended ~ 6 kpc and 1 
kpc, respectively, which are 6-12 times the sizes 
of the respective CO distributions listed in Table 
4. This issue will be addressed further in a follow- 
up paper containing CO interferometric maps of 
the sample (Evans et al. 2005a). 

5. DISCUSSION 

The present CO survey has increased the num- 
ber of radio galaxies observed in Paper 1 by a 
factor of four. Three new CO detections have 
been made, yielding a 25% detection efficiency 
for the sample with the 12m telescope. The CO 



detections translate into molecular gas masses of 
0.4 — 7 x 10 9 M , with 3<r upper limits for non- 
detections in the range of < 4 x 10 8 M Q for a 
z ~ 0.02 galaxy to < 3.3 x 10 9 M for a z ~ 0.1 
galaxy. Thus, while only a quarter of the sam- 
ple was detected, the upper limits of the more 
distant members of the sample do not rule out 
molecular gas mass as high as that of the Milky 
Way (~3x 10 9 M ). 

There are two notable features in the CO spec- 
tra shown in Figure 1. The first is the broad 
linewidth - the average Aw FW hm of the present 
sample of radio galaxies is ~ 500 ± 130 km s _1 , 
which exceeds the average value for both PG 
QSOs (260±160 km s' 1 : Evans et al. 2005b) 
and ultraluminous infrared galaxies (300±90 km 
s _1 : Solomon et al. 1997) detected in CO to 
date. This result likely indicates that the host 
galaxies of molecular gas-rich radio galaxies are 
extremely massive. The second notable feature 
is the double-peaked CO emission line profiles in 
many of the spectra. Interferometric CO observa- 
tions of two of these galaxies are in the literature 
and show this feature to be present for (at least) 
two different reasons. In the case of 3C 293, the 
"double-peak" appearance is caused by a CO ab- 
sorption feature at the systemic velocity of the 
galaxy (Evans et al. 1999a). In the case of 3C 31, 
the molecular gas is distributed in a ring (Okuda 
et al. 2005). 

Figure 4 shows optical images of 7 of the nine 
radio galaxies detected in CO; two of the nine 
radio galaxies (B2 1318+34 and PKS 1345+12) 
are well-studied systems that have been imaged 
by several groups (e.g. Sanders et al. 1988b; 
Bushouse & Stanford 1992; Surace et al. 1998; 
Evans et al. 1999b; Scoville et al. 2000; Kim, 
Veilleux, & Sanders 2003), and thus have not 
been reobserved by us. Most of these molecu- 
lar gas-rich galaxies have nearby companions (3C 
31; PKS 1345+12) and/or show clear evidence 
of being involved in a recent merger (3C 84; 3C 
120; PKS 0502-10; 3C 293; PKS 1345+12; B2 
1318+34). The one exception is B2 0722+30, an 
edge-on spiral galaxy. While the images in Fig- 
ure 4 support the possibility that interaction and 
merger events may be a primary trigger for radio- 
loud AGN activity, it is notable that many low- 
redshift radio galaxies are elliptical galaxies (e.g., 
Sandage 1966; Martel et al. 1999). Thus, unless 
the elliptical galaxies are advanced mergers, other 
mechanisms must play an important role in ignit- 
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ing the AGN activity observed in these galaxies. 

Given the results of the present survey and re- 
cent CO surveys of low-redshift PG QSOs (Evans 
et al. 2001; Scoville 2003), there is sufficient data 
to begin a comparison of the molecular gas prop- 
erties and the global energy output (measured at 
infrared, radio, and X-ray wavelengths) of clas- 
sical UV and radio-selected AGN with infrared- 
selected galaxies and elliptical galaxies. Before 
exploring such comparisons, it is worth first con- 
sidering two issues which have direct relevance to 
the interpretation of the data. The first issue is 
the infrared bias in the sample. While the sam- 
ple was selected based on the 60/xm and 100/xm 
flux densities, it is important to realize that (a) 
the flux density limits are extremely low, and 
thus the infrared luminosity range of the sample 
spans several decades (i.e., Lir ~ 10 9-12 L ), and 
that (b) 55%, 33%, and 24% of z £ 0.1 3C, B2 
and PKS galaxies, respectively, were detected at 
60/xm by IRAS. The fraction for PG QSOs is even 
higher, with 69% of PG QSOs out to z ~ 0.17 
(the redshift limit of the QSO CO survey) hav- 
ing 60/zm detections. Thus, infrared excesses in 
AGN host galaxies from radio-selected (and UV- 
selected) surveys are quite common. 

The second issue of relevance to the interpreta- 
tion of the data is the nature of the radio emis- 
sion from the galaxies in the sample; i.e., whether 
the radio emission from all of the galaxies in the 
sample is due to AGN jet emission, or whether 
some of the galaxies are actually nearby star- 
burst galaxies luminous enough in the radio to 
have been detected by the 3C, B2 or PKS sur- 
veys. To address this issue, IRAS and 1.4 GHz 
data have been compiled for the radio galax- 
ies, along with complementary data sets for in- 
frared galaxies with both warm and cool 25-to- 
60/im colors (note that classical AGN typically 
have warm colors, /^m/ieo/um > 0.2), PG QSOs, 
and elliptical galaxies. Figure 5 shows the re- 
sult of the compilation; plotted is the 1.4 GHz-to- 
infrared luminosity ratio, L 1A q Hz , versus the in- 
frared luminosity, L IR . The "radio- infrared" cor- 
relation (see Condon 1992) is clearly visible in this 
Figure; the radio-to-infrared ratio for infrared- 
selected galaxies is relatively constant with in- 
creasing infrared luminosity. This has led to spec- 
ulation that the same population of massive stars 
responsible for heating the dust is also responsi- 
ble for the radio emission (through, for example, 
the production of supernovae). Lonsdale, Smith, 



& Lonsdale (1995) showed that QSOs follow the 
same bolometric luminosity-radio correlation as 
infrared galaxies. Similarly, the QSOs follow the 
infrared-radio correlation in Figure 5, indicating, 
perhaps, that a significant fraction of their ra- 
dio and infrared emission is associated with mas- 
sive stars, or that AGN-related processes in radio- 
quiet sources yield the proper Li.4gh z /-^ir ratio 
to place QSOs on the relation. In contrast, the 
data points for all but one of the radio galaxies 
are above this correlation; their excess radio emis- 
sion is associated with AGN jets. Only the galaxy 
B2 1318+34 (aka IC 883), which is classified op- 
tically as a LINER galaxy, has a Li.4ghz/ ^ir ra- 
tio consistent with the infrared-radio correlation. 
The LINER classification is commonly an indica- 
tion of ionization of an AGN with a shallow power 
law, however, the spatial coincidence between the 
radio (Condon et al. 1991) and CO (Downes & 
Solomon 1998) emission in the disk of IC 883 is 
evidence in favor of the radio emission being as- 
sociated with a starburst population. Given this, 
B2 1318+34 will be denoted by a symbol different 
from the other radio galaxies in Figures 6 and 8. 

5.1. Global CO -Infrared Properties 

Figure 6 shows two common plots of infrared 
and CO luminosity for extragalactic sources. The 
first plot is the CO luminosity versus the infrared 
luminosity of the radio galaxies and PG QSOs, 
plotted with infrared galaxies and elliptical galax- 
ies with CO detections. In addition, arrows have 
been plotted to represent 3a L' co upper limits for 
the radio galaxies with CO nondetections in the 
present survey. There are two noteworthy fea- 
tures in Figure 6a. First, the PG QSO population 
appear to have low CO luminosities for their in- 
frared luminosity relative to infrared galaxies (see 
also Evans et al. 2001). Second, by contrast, the 
radio galaxies are well distributed throughout the 
area of the plot occupied by infrared and ellipti- 
cal galaxies, with the three newly detected radio 
galaxies (3C 31, PKS 0502-10, and 3C 293) hav- 
ing high CO luminosities for their infrared lumi- 
nosity. These features are better represented in 
Figure 6b, which is a plot of the L IR /L co ver- 
sus the infrared luminosity. The ratio Lm/ L' co 
is commonly referred to as the star formation 
efficiency, i.e., it is the measure of energy from 
massive stars heating the dust normalized by the 
amount of molecular gas available to form new 
stars. Thus, for highly efficient starbursts (i.e., 
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high Lir/Lq ), a large number of massive stars 
are produced from the available gas. The QSOs 
have high Lir/L co , indicating that either they 
are producing massive stars with extremely high 
efficiency, or that there is significant heating of the 
dust by the AGN (Evans et al. 2001). Most of the 
radio galaxies have Li^/L' co similar to those of gi- 
ant molecular clouds (GMCs: Sanders & Mirabel 
1996) in the Galaxy. Exceptions to this include 
PKS 1345+12, which has a Li R /Z/ co similar to 
the PG QSOs, and the three newly detected ra- 
dio galaxies, which have L IR /L co (< 30) similar 
to the global value for local spiral galaxies (see 
also Evans et al. 1999a). 

It thus appears that while some radio galax- 
ies have significant amounts of molecular gas, the 
rate of massive star formation is extremely low. 
Indeed, if the assumption is made that the in- 
frared luminosity in 3C 31, PKS 0502-10, and 3C 
293 is due entirely to star formation, then the star 
formation rates are roughly a few solar masses per 
year. Note also that the dust temperatures of 3C 
31, PKS 0502-10, and 3C 293 are lower (30-45 K) 
than that of the previously detected radio galaxies 
(55-90 K), also consistent with the idea that cur- 
rent massive star formation rates in these three 
radio galaxies is low. 

5.2. Global X-ray Properties 

Given the amount of molecular gas in the radio 
galaxies detected in CO to date, and the diversity 
in the apparent levels of active star formation (as 
traced by Li R and Lir/L' co ), it is worth examin- 
ing their hard (2-10 keV) X-ray properties to de- 
termine (a) whether these radio galaxies have X- 
ray emission indicative of AGN, and (b) whether 
any of the radio galaxies have sufficient column 
densities of molecular gas to be significantly X- 
ray absorbed. The Lx-to-Lpi R ratio provides a 
straight-forward way of addressing these issues - 
the ratio is constant in starburst galaxies, indicat- 
ing that the starburst responsible for heating the 
dust in HII region galaxies produces a constant 
fraction of massive X-ray binary pairs (Grimm et 
al. 2003; Ranalli et al. 2003). 'Excesses' in L x 
I Lfir arise from the contribution of AGN to the 
X-ray emission, making the ratio a powerful AGN 
vs. starburst diagnostic. However, in objects that 
contain large columns of gas and dust, the X-ray 
emission can be suppressed, lowering the observed 

Lx/ -^FIR- 

For the purposes of the present discussion, the 



available radio galaxy data of our sample have 
been compiled from the literature and added to 
the Lx vs. Lfir plot of Trentham et al. (2005), 
which is a compilation of data of local LINERs, 
Seyferts, Hll-region galaxies, infrared galaxies, 
and PG QSOs. The result is plotted in Figure 7; 
Figure 7a shows the data plotted in the manner 
typically done by other authors, i.e., Lx vs. Lfir, 
and Figure 7b shows the data plotted in terms 
of the Lx/£fir versus Lfir- A least-squares fit 
has been applied to the HII region galaxy data. 
Several of the infrared galaxies, which are rich in 
molecular gas, are below the fit and thus Comp- 
ton thick. In terms of the radio galaxies data, the 
Figures show that, (a) the AGN dominates the X- 
ray emission from the radio galaxies in a manner 
consistent with PG QSOs and lower luminosity 
Seyferts and LINERs, and (b ) the X-ray emission 
from the AGN does not seem to be significantly 
suppressed, and thus these radio galaxies do not 
reside in Compton-thick environments. Indeed, 
several of the radio galaxies have measured col- 
umn densities in the literature - the derived col- 
umn densities from these X-ray observations are 
~ 10 21 - 22 cm.- 2 (Leighly et al. 1997; Sambruna 
et al. 1999; and Odea et al. 2000). 

5.3. CO Luminosity and L IR /L' co versus Radio 

Power 

Another issue relevant to the present sample of 
galaxies is whether the molecular gas properties 
are correlated with the radio output of the radio 
galaxies. I.e., if molecular gas fuels AGN activ- 
ity, there may be a correlation between L' co and 
the energy emitted by these AGN at radio wave- 
lengths. Figure 8a is a plot of L' co versus Li.4gh z 
of the sample of radio galaxies; comparative sam- 
ples of infrared galaxies, PG QSOs, and ellipti- 
cals are also plotted. The infrared galaxy and PG 
QSO data are correlated, with higher L' co corre- 
sponding to higher Li.4gh z - Such a correlation can 
be understood, for example, if a significant frac- 
tion of the 1.4 GHz emission is generated via star 
formation - the higher the molecular gas content, 
the higher the number of supernovae generating 
radio emission. The radio galaxies occupy a differ- 
ent location on this plot, having high Li.4gh z f° r 
their L' co relative to infrared galaxies and QSOs. 
It is also clear from the radio galaxy data that 
no obvious correlation exists between Li.4gh z and 
L' co . This result may simply illustrate the fact 
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that, if a correlation between the AGN radio out- 
put and molecular gas mass exists, it may only 
exist for the component of the molecular gas in 
the nuclear region (near the AGN), and not the 
global molecular gas mass measured by these ob- 
servations. 

Figure 8b is a plot of Lir./ L' co versus Li.4gh z - 
In terms of the distribution of infrared galaxies 
and QSO data, the plot is reminiscent of Figure 
6b - L\Yi/L' co increases with increasing Li.4gh z , 
and the QSOs and warm infrared galaxies reside 
primarily at high Li K /L' co and L 1A q Hz relative 
to cool infrared galaxies. The radio galaxies oc- 
cupy a different location on the plot, having high 
Liaghz for their Lm/L' co , but no obvious corre- 
lation among the radio galaxy data exists. Thus, 
regardless of whether stars or AGN in the host 
galaxy of these radio galaxies are heating the dust, 
the fueling efficiency factor Li&/ L' co is not corre- 
lated with the total radio output. 

5.4. High Mass Star Formation 

All but one of the radio galaxies detected as 
part of this survey have molecular gas masses in 
excess of 1 x 10 9 M , indicating that there are 
large fuel reservoirs for star formation and pos- 
sibly AGN activity. The presence of molecular 
gas, in itself, is not evidence of active star for- 
mation. High mass star formation in the Milky 
Way occurs in the cores of giant molecular clouds 
(GMCs) where the molecular gas densities exceed 
~ 10 5 cm -3 . 

Preliminary observations have been done of 
PKS 0502-10 and B2 0722+30 to provide addi- 
tional information about the state of the molecu- 
lar gas. The first observations are of the CO (2 — > 
1) emission line; the -^co(2-»i)/-^co(i-»o) r& tios of 
PKS 0502-10 and B2 0722+30 are 0.6 and 0.7, re- 
spectively. Such low values of -^co(2-»i)/^co(i^o) 
are suggestive of subthermally excited CO, i.e., 
low molecular gas temperatures and densities. 

The second set of preliminary observations de- 
signed to probe the state of the gas in these radio 
galaxies are of the high density tracer HCN. HCN 
is seen in abundance when the molecular gas den- 
sities are on the order of 10 5 cm -3 , and thus is 
found in the cores of GMCs where massive star 
formation is commencing. Gao & Solomon (2004) 
have presented an HCN survey of infrared galax- 
ies, and have calculated an average Lir/L^ cn ra- 
tio of ~ 970 ± 600. The HCN(1 -> 0) obser- 
vations of two radio galaxies yielded only upper 



limits, translating into Lir/L^ cn lower limits of ~ 
550 and 3700 for PKS 0502-10 and B2 0722+30, 
respectively. Thus, while the Lir/L^ cn of PKS 
0502-10 does not provide a useful constraint, the 
high Li R /L' RCN of B2 0722-30 is suggestive of ei- 
ther highly efficient massive star formation or sig- 
nificant heating of dust by the AGN. 

5.5. Are Fan aroff- Riley II Radio Galaxies 
Relatively Molecular Gas Poor? 

One of the goals of this survey was to obtain 
detections of CO emission from radio galaxies of 
all morphological types. This survey yielded a 
CO detection in only one Fanaroff- Riley II (edge- 
brightened radio morphology: Fanaroff & Riley 
1974) radio galaxy, 3C 293. However, in terms 
of its radio properties, 3C 293 is not a represen- 
tative F-R II. The galaxy has a low radio power 
(£i.4GHz ~ 2 x 10 25 W m" 2 Hz" 1 ) relative to the 
bulk of F-R II radio galaxies. In addition, 3C 293 
has a complex steep spectrum core which shows 
evidence of a previous radio outburst (Akujor et 
al. 1996). 

The low CO detection rate of F-R II galaxies (~ 
8%) relative to F-R I and compact radio source 
galaxies (~ 35%) may be a simple consequence 
of larger distances of F-R II galaxies, on average, 
relative to lower radio power radio galaxies in the 
sample. This explanation is unlikely; Figure 6a 
shows the L' co of the detected radio galaxies and 
the L' co upper limits of the undetected galaxies. 
The upper limits of all of the non-detected F-R 
II galaxies are below 2.5 x 10 9 K km s _1 pc 2 , 
and three out of seven (excluding B2 1318+34) of 
the CO-detected F-R I and compact radio galax- 
ies are above this line. Thus, though the F-R II 
galaxies comprise 36% of the sample, there are 
no additional F-R II galaxies in the sample with 
molecular gas masses comparable to either 3C 
293, or the most gas-rich F-R I and radio compact 
galaxies (PKS 0502-10, TXS 1506+345, and PKS 
1345+12) in the sample. This result, of course, 
relies on the assumption that F-R II radio galax- 
ies have CO emission line widths of ~ 500 km 
s^ 1 ; if the line widths are sufficiently broad, then 
the limited velocity bandwidth (i.e., ~ 1400 km 
s _1 at 3mm) available with single-dish telescopes 
such as the 12m Telescope would make the CO 
detections difficult. 

The present CO survey has focussed on a 
sample of IRAS-detected radio galaxies using a 



10 



CO IN RADIO GALAXIES 



fairly moderate size single-dish millimeter tele- 
scope. The results presented here can be im- 
proved upon in the following ways: (i) CO sur- 
veys of larger samples of ISO and Spitzer Space 
Telescope-detected radio galaxies using larger 
single-dish millimeter-wave telescopes or interfer- 
ometers, (ii) Volume-limited CO surveys of ra- 
dio galaxies (and QSOs) similar to those initi- 
ated by Lim et al. (1999) and Scoville et al. 
(2003). Such surveys will provide a more com- 
plete view of the molecular gas properties of ra- 
dio (and UV-selected) AGN. (Hi) Interferometric 
observations of the CO-detected radio galaxies to 
measure the extent, surface density, and kinemat- 
ics of the molecular gas relative to the stellar pop- 
ulation and the ionized/neutral gas (Inoue et al. 
1996; Evans et al. 1999a,b; Okuda et al. 2005; 
Evans et al. 2005a). Finally, continued investi- 
gations of star formation in radio galaxies using 
diagnostics at multiple wavelengths will provide a 
more complete picture of the putative starburst- 
AGN connection. Recent optical spectroscopy of 
three of the galaxies in the sample (3C 293, 3C 
305, and PKS 1345+12: Tadhunter et al. 2005) 
have provided encouraging results. 

6. SUMMARY 

The survey presented here is a follow-up to the 
Mazzarella et al. (1993) study of IRAS-detected 
radio galaxies, which yielded new CO detections 
in three radio galaxies. The survey has increased 
the number of low-redshift radio galaxies observed 
in the prior survey by a factor of four. The fol- 
lowing conclusions have been reached: 

(1) Three additional radio galaxies (3C 31, PKS 
0502-10, 3C 293) have been detected in CO(l -> 
0), yielding a total of 9 detections out of 33 radio 
galaxies. The CO detections translate into molec- 
ular gas masses of 0.4 — 7 x 10 9 M , with 3a upper 
limits for nondetections in the range of < 4 x 10 s 
M for a z ~ 0.02 galaxy to < 3.3 x 10 9 M for 
a z ~ 0.1 galaxy (a = 1.5 M [K km s _1 pc 2 ] -1 . 
Thus, the upper limits of the more distant mem- 
bers of the sample do not rule out molecular gas 
mass as high as that of the Milky Way (~ 3 x 10 9 
M ). 

(2) Many of the CO lines have double-peaked 
profiles, which is indicative of either CO absorp- 
tion towards the AGN, or CO distributed in a 
molecular ring. 

(3) The CO line widths are extremely broad. 
The average velocity dispersion of the CO- 



detected radio galaxies is Aw FW hm = 500 ± 130 
km s" 1 , which exceeds the average line width of 
both ultraluminous infrared galaxies (= 300 ± 90 
km s" 1 ) and PG QSOs (= 260 ± 160 km s" 1 ). 
This is an indication that the CO-luminous radio 
galaxies are massive. 

(4) With the exception of the edge-on spi- 
ral galaxy B2 0722+30, optical images show the 
molecular gas-rich radio galaxies to have close 
companions and/or tidal features consistent with 
those of infrared- luminous interacting and merg- 
ing galaxy systems. 

(5) Most of the radio galaxies have L IR /L' co 
similar to GMCs. Exceptions to this include (i) 
PKS 1345+12, which has a L IR /L' co (~ 370) com- 
parable to the PG QSOs, likely an indication 
of highly efficient massive star production or of 
significant dust heating by the AGN, and (ii) 
3C 31, PKS 0502-10, and 3C 293, which have 
low Lir/L' co (< 30) consistent with the global 
L IR /Z/ co of spiral galaxies. 

(6) All of the radio galaxies have measured val- 
ues of Z/x/ -Zvfir significantly higher than the aver- 
age for HH-region like galaxies. Thus, the AGN is 
the major contributor to the X-ray emission, and 
the AGN do not reside in Compton-thick environ- 
ments. 

(7) No correlation exists between L' co or 
L m / L' co and Li.4gh z of the radio galaxies. 

(8) Only one F-R II radio galaxy (3C 293) was 
detected, yielding an 8% CO detection rate for 
F-R II galaxies. By comparison, ~ 35% of the 
F-R I and compact radio source galaxies were de- 
tected. The L' co upper limits of the undetected F- 
R II galaxies are sufficiently small to support the 
conclusion that either low-redshift F-R II galax- 
ies are gas-poor compared with F-R I and radio 
compact galaxies, or that the CO line widths of 
the undetected F-R II galaxies are broad enough 
to make their detection difficult with a single-dish 
telescope. 
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Figure Captions 



Figure 1. (a) NRAO 12m and IRAM 30m CO(l— >0) spectrum of 31 of the 33 radio galaxies in the 
survey. The intensity scale is in units of main beam brightness temperature. A linear baseline has been 
subtracted from each spectrum; the baseline subtraction is performed outside of the velocity range of 
emission lines when they are present. With the exception of PKS 0502-10, the zero velocity corresponds 
to the redshift of the radio galaxy listed in Table 1. The zero velocity redshift of PKS 0502-10 is 0.0396. 
(b) IRAM 30m CO(l-» 0) spectrum of PKS 1345+12. The spectrum confirm the general features of 
the 12m spectrum of Figure la, but with significantly higher signal-to-noise. 

Figure 2. JCMT 15m and IRAM 30m telescope CO(2 -> 1) spectra of PKS 0502-10, B2 0722+30, and 
Cygnus A. The intensity scale is in units of main beam brightness temperature (IRAM) and antenna 
temperature corrected for aperture losses (JCMT). A linear baseline has been subtracted from each 
spectrum; the baseline subtraction is performed outside of the velocity range of emission lines when 
they are present. Zero velocity corresponds to the redshift of the radio galaxy listed in Table 3. 
Figure 3. IRAM 30m HCN(1 -> 0) spectra of PKS 0502-10 and B2 0722+30. The intensity scale is in 
units of main beam brightness temperature. A linear baseline has been subtracted from each spectrum; 
the baseline subtraction is performed outside of the velocity range where the HCN emission line is 
expected to be. Zero velocity corresponds to the redshift of the radio galaxy listed in Table 3. 
Figure 4. False-color U , B, and /-band images of seven of the radio galaxies detected in CO(l — > 0). 
The U, B, and /-band data are displayed as blue, green, and red, respectively. North is up, and east is 
to the left. 

Figure 5. The ratio of 1.4 GHz to infrared luminosity versus the infrared luminosity of infrared luminous 
radio galaxies, a sample of infrared galaxies, PG QSOs and elliptical galaxies detected in CO. References 
for the infrared data are as follows: infrared galaxies, Mazzarella et al. (1993); PG QSOs, Sanders et 
al. (1989) and Haas et al. (2003); elliptical galaxies (Moshir et al. 1990). References for the 1.4 GHz 
data are as follows: infrared galaxies, Condon et al. (1990, 1991, 1992); elliptical galaxies (White & 
Becker 1992; Becker et al. 1995; Condon et al. 1998); radio galaxies (White & Becker 1992; Becker 
et al. 1995); PG QSOs (White & Becker 1992; Barvainis, Lonsdale, & Antonucci 1996; Condon et al. 
1998). 

Figure 6. (a) CO vs. infrared luminosity plot of infrared luminous radio galaxies, a sample of infrared 
galaxies, PG QSOs and elliptical galaxies detected in CO. References for CO data are as follows: infrared 
galaxies, Mazzarella et al. (1993); PG QSOs, Sanders et al. (1988c), Barvainis et al. (1989), Evans et 
al. (2001), and Scoville et al. (2003); elliptical galaxies: Wiklind & Combes (1997). (b) the ratio of 
infrared to CO luminosity vs. infrared luminosity for the same sample of galaxies as in 6(a). In both 
Figures, the galaxy B2 1318+34 is denoted by an unfilled triangle. 

Figure 7. (a) Hard (2-10 keV) X-ray luminosity vs. far-infrared luminosity plot of the infrared luminous 
radio galaxies and a comparison sample of infrared galaxies, PG QSOs, and less luminous galaxies 
classified as HII region-like galaxies, Seyferts, and LINERs. The far-infrared luminosities have been 
calculated using data in NED. The radio galaxies in the Figure are NGC 315, 3C 84, 3C 111, 3C 
120, PKS 0634-20, Hydra A, M87, PKS 1345+12, 3C 293, 3C 321, NGC 6251, 3C 390.3, and 3C 405. 
References for the X-ray data are as follows: HII, Seyferts, and LINERs, Terashima et al. (2002) and 
Grimm et al. (2003); infrared luminous galaxies, Iwasawa & Comastri (1998), Boiler (2002), Dalla Ceca 
(2002), Franceschini et al. (2003), and Sanders et al. (2005); radio galaxies, Churazov et al. (2003); Di 
Matteo et al. (2003); Leighly et al. (1997); McNamara et al. (2000); Odea et al. (2000); and Sambruna 
et al. (1999). The X-ray luminosity of 3C 293 (Lx ~ 3x 10 42 erg s _1 ) was calculated by us from archival 
ASCA data, (b) Ratio of the X-ray to far-infrared luminosity versus the far-infrared luminosity for the 
same sample of galaxies as in Figure 7(a). In both (a) and (b), radio galaxies and PG QSOs with CO 
detections are encircled. Adapted from Trentham et al. (2005). 

Figure 8. (a) The CO luminosity vs. the radio power of the infrared luminous radio galaxies and a 
comparison sample of infrared galaxies, PG QSOs, and elliptical galaxies, (b) Lfr/ L' co vs. the radio 
power for the same sample of galaxies plotted in 8a. In both Figures, the galaxy B2 1318+34 is denoted 



by an unfilled triangle. 



Table 1 
Source List 





Coordinates (J2000.0) 






./60/xm 


/lOOpm 


ilR a 


iFIH a 


Pl.5GHz 


Source 


RA 


Dec 


z 


FR 


(mjy) 


(mjy) 


( L 0) 


( L <s>) 


(W m~ 2 Hz" 1 ) 


NGC 315 


00h57m48.85s 


+30d21m08.63s 


0.0165 


I 


368 


460 


7.2 xlO 9 


2.3xl0 9 


1.6X10 24 


3C 31 


01h07m24.98s 


+32d24m44.77s 


0.0167 


I 


435 


1675 


7.4xl0 9 


4.7X10 9 


2.5xl0 24 


4C 31.04 


01hl9m35.00s 


+32dl0m50.03s 


0.0590 


C 


150 


524 


3.7xl0 1() 


1.9X10 10 


1.7xl0 25 


NGC 741 


01h56m21.05s 


+05d37m44.17s 


0.0185 


C 


214 


777 


7.0xl0 9 


2.7X10 9 


6.0 xlO 23 


3C 84 


03hl9m48.09s 


+41d30m42.47s 


0.0176 


I 


7427 


8267 


1.6X10 11 


5.1X10 10 


7.3xl0 24 


3C 88 


03h27m54.17s 


+02d33m41.85s 


0.0302 


11 


180 


816 


1.7xl0 1() 


7.0X10 9 


8.5xl0 24 


3C 111 


04hl8m21.07s 


+38d01m32.59s 


0.0485 


11 


321 


<2250 


9.0xl0 10 


2.8X10 10 


6.5xl0 25 


3C 120 


04h33mll.09s 


+05d21ml6.03s 


0.0331 


I 


1383 


1937 


1.2X10 11 


3.6X10 10 


9.0xl0 24 


PKS 0502-10 


05h04m52.99s 


-10dl4m52.18s 


0.0410 


? 


688 


1370 


6.1xl0 10 


3.0X10 10 


4.6xl0 24 


PKS 0634-20 


06h36m32.24s 


-20d34m52.92s 


0.0544 


II 


484 


1150 


9.5xl0 10 


4.3X10 10 


3.3xl0 25 


B2 0648+27 


06h42m02.82s 


+27d28m21.69s 


0.0414 


c 


2633 


1529 


2.2X10 11 


8.7X10 10 


4.9xl0 23 


B2 0722+30 


07h25m37.35s 


+29d57ml4.71s 


0.0188 


I 


3108 


4999 


5.2xl0 1() 


2.8X10 10 


8.2xl0 22 


4C 29.30 


08h40m02.44s 


+29d49m02.87s 


0.0650 


I 


472 


595 


l.lxlO 11 


4.7X10 10 


3.3xl0 23 


Hydra A 


09hl8m05.69s 


-12d05m45.45s 


0.0538 


I 


155 


416 


2.8xl0 10 


1.4xl0 10 


2.5xl0 26 


3C272.1 


12h25m03.78s 


+12d53ml3.1s 


0.0035 


I 


556 


1024 


5.6xl0 8 


1.8x10 s 


l.OxlO 23 


M 87 


12h30m49.35s 


+12d23m28.06s 


0.0047 


I 


546 


559 


l.lxlO 9 


2.6x10 s 


7.2xl0 24 


B2 1318+34 


13h20m35.20s 


+34d08ml3.56s 


0.0233 


I 


16070 


20682 


3.4x10" 


2.0X10 11 


8.9xl0 22 


3C 285 


13h21ml7.78s 


+42d35ml4.23s 


0.0789 


II 


253 


637 


1.8x10" 


5.0xl0 10 


2.6xl0 25 


PKS 1345+12 


13h47m33.38s 


+12dl7m24.18s 


0.1224 


c 


2098 


1738 


1.5xl0 12 


6.7X10 11 


1.5xl0 26 


3C 293 


13h52ml7.86s 


+31d26m46.67s 


0.0450 


II 


233 


621 


2.9xl0 10 


1.5xl0 10 


1.7 xlO 25 


OQ 208 


14h07m00.39s 


+28d27ml4.73s 


0.0768 


c 


753 


1029 


4.1X10 11 


l.lxlO 11 


8.3xl0 24 


3C 305 


14h49m21.46s 


+63dl6ml6.14s 


0.0410 


I 


298 


450 


2.4xl0 ll) 


1.2xl0 10 


9.4X10 24 


TXS 1506+345 


15h08m05.64s 


+34d23m22.75s 


0.0449 


I 


2647 


3956 


2.3x10" 


1.3X10 11 


4.5 xlO 23 


PKS B1518+045 


15h21m22.61s 


+04d20m29.21s 


0.0513 


? 


<120 


968 


4.6xl0 1() 


1.8xl0 10 


2.1xl0 25 


3C 321 


15h31m43.45s 


+24d04ml8.70s 


0.0960 


II 


1067 


961 


5.0X10 11 


2.1X10 11 


6.2xl0 25 


3C 327 


16h02m27.39s 


+01d57m55.70s 


0.1048 


II 


670 


371 


4.1X10 11 


1.4X10 11 


1.9xl0 26 


NGC 6251 


16h32m31.97s 


+82d32ml6.46s 


0.0249 


I 


188 


600 


5.7xl0 9 


2.7xl0 9 


2.6xl0 24 


B2 1707+34 


17h09m38.37s 


+34d25m53.76s 


0.0808 


II 


464 


841 


2.0X10 11 


8.2xl0 10 


8.4 xlO 24 


3C 390.3 


18h42m08.92s 


+79d46ml7.27s 


0.0561 


11 


250 


340 


1.3X10 11 


1.9xl0 10 


7. 2 xlO 25 


3C 402 


19h41m46.00s 


+50d35m44.86s 


0.0254 


II 


257 


1052 


l.lxlO 10 


6.6 xlO 9 


3.1 xlO 24 


3C 403 


19h52ml5.36s 


+02d30m28.14s 


0.0590 


II 


441 


<500 


1.5X10 11 


3.0xl0 10 


3.9 xlO 25 


3C 405 


19h59m28.51s 


+40d44m01.70s 


0.0556 


11 


2847 


<1800 


4.0X10 11 


1.6X10 11 


9.4xl0 27 


3C 433 


21h23m44.53s 


+25d04m09.99s 


0.1016 


I 


299 


<1200 


2.9X10 11 


8.8xl0 10 


2.3xl0 24 



a Both Lm and ^FIR are calculated using the prescription in Sanders & Mirabel (1996). In cases where one or more IRAS flux 
densities have only upper limit values, LrR and Lfir are calculated by setting the flux density value equal to the upper limit value and 
by setting the flux density value equal to zero, then taking the average ol these two determinations. 

References. — Fanaroff-Rilcy classications: Golombek et al. 1988; Baum, Heckman, & van 

Brcugel 1992; http://www.jb.man.ac.uk/atlas/. Flux Densities: Golombek et al. 1988; Rcdshifts: Golombek ct al. 1988; Mirabel 
et al. 1989; Mazzarella et al. 1993. 



Table 2 
Journal of Observations 



Date 



Time 



Source 



Site Transition (mmm-yy) (hr) 



T 

J- sys 

(K) 



Oct-97 


3.0 


420 


Jan-96 


17.4 


272 


Jun-97 


18.4 


333 


Oct-97 


4.6 


273 


Jan-96 


29.6 


199 


Jun-97 


12.4 


371 


Jun-97 


4.8 


318 


Oct-97 


7.2 


248 


Sep-95 


10.2 


260 


Oct-97 


19.6 


259 


Sep-95 


14.8 


369 


Jan-96 


9.4 


256 


Sep-95 


4.0 


411 


Dec-03 


1.8 


120 


Jun-97 


5.2 


372 


Oct-97 


12.6 


256 


Sep-95 


24.6 


319 


Jan-96 


2.6 


237 


Oct-97 


6.6 


299 


Sep-95 


1.8 


469 


Dec-03 


1.8 


83 


Jan-96 


12.2 


220 


Jun-97 


30.8 


232 


Oct-97 


8.2 


196 


Jun-97 


5.2 


288 


Oct-97 


9.4 


224 


Dec-92 


21.9 


368 


Oct-97 


1.8 


256 


Jan-96 


49.6 


210 


Jun-97 


30.1 


256 


Jun-00 


4.3 


200 


Jan-96 


17.6 


240 


SeD-95 


12.8 


261 


Oct-97 


16.0 


261 


Oct-97 


8.4 


211 


Jun-97 


9.4 


269 


Sep-95 


23.6 


276 


Jun-97 


16.4 


249 


Oct-97 


5.6 


321 


Jun-00 


16.2 


205 


Jun-97 


13.0 


252 


Oct-97 


12.0 


208 


Jun-97 


16.4 


356 


Oct-97 


5.2 


464 


Oct-97 


6.0 


248 



NGC 315 
3C 31 



4C 31.04 
NGC 741 
3C 84 
3C 88 
3C 111 
3C 120 
PKS 0502-10 



PKS 0634-20 
B2 0648+27 
B2 0722+30 

4C 29.30 

Hydra A 
M 87 

B2 1318+34 
3C 285 

PKS 1345+12 

3C 293 
OQ 208 
3C 305 

TXS 1506+345 
PKS B1518+045 
3C 321 

3C 327 

B2 1707+34 

3C 390.3 

3C 402 



12m 
12m 
12m 
12m 
12m 
12m 
12m 
12m 
12m 
12m 
12m 
12m 
15m 
30m 
12m 
12m 
12m 
12m 
12m 
15m 
30m 
12m 
12m 
12m 
12m 
12m 
12m 
12m 
12m 
12m 
12m 
12m 
12m 
12m 
12m 
12m 
12m 
12m 
12m 
12m 
12m 
12m 
12m 
12m 
12m 



CO( 
CO( 
CO( 
CO( 
CO( 
CO( 
CO( 
CO( 
CO( 
CO( 
CO( 
CO( 
CO( 
HCN(1- 



CO( 
CO( 
CO( 
CO( 
CO( 
CO( 
HCN(1- 



co 

CO 
CO 
CO 
CO 
CO 
CO 
CO 
CO 
CO 
CO 
CO 
CO 
CO 
CO 
CO 
CO 
CO 
CO 
CO 
CO 
CO 
CO 
CO 



0) 



0) 
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Table 2 — Continued 



Date Time T sys 

Source Site Transition (mmm-yy) (hr) (K) 

3C 403 12m CO(l-> 0) Jun-97 11.4 254 

3C 405 30m CO(l-> 0) Nov-95 5.9 186 

30m CO(2^ 1) Nov-95 3.3 321 

3C 433 12m CO(l-> 0) Jan-96 11.6 200 

12m CO(l-^O) Oct-97 13.0 179 



Table 3 

Molecular Gas Emission Line Properties 



Source 


n a 
L 


line 


OltC 


ZCO 


AllpwHM 


T 1 A o. b 


•^lincAl) 


Li,. 


1V1 \ri-2 ) 




(MpcJ 










(km s ) 


(rv km s ) 


(Jy km s~ 1 ) 


1 9 

(K km s~ pc ) 




NGC 315 


66 


CO(l- 


.0) 


12m 






<0.91 


<23.0 


<2.4xl0 8 


<3. 6x10 s 


3C 31 


67 


CO(1->0) 


12m 


0.0169 


450 


1.1+0.08 


27+2 


2.9xl0 8 


4.3xl0 8 


4C 31.04 


239 


CO(l- 


.0) 


12m 






<0.21 


<5.2 


<6.7xl0 8 


<1.0xl0 9 


NGC 741 


74 


CO(1^0) 


12m 






<0.39 


<9.8 


<1.3xl0 8 


<1. 9x10 s 


3C 84 


71 


CO(l- 


>0) 


12m 


0.0176 


200 


4.1+0.05 


104+1 


1.2xl0 9 


1.8xl0 9 


3C 88 


122 


CO(l- 


>0) 


12m 






<0.41 


<10.4 


<3.6xl0 8 


<5. 4x10 s 


3C 111 


196 


CO(l-vO) 


12m 






<0.34 


<8.5 


<7.5xl0 8 


<l.lxl0 9 


3C 120 


133 


CO(l-vO) 


12m 


0.0331 


550 


1.1+0.1 


28+2 


1.2X10 9 


1.7xl0 9 


PKS 0502-10 


160 


CO(l-vO) 


12m 


0.0398 


550 


1.8+0.09 


45+2 


2.7xl0 9 


4.0xl0 9 






CO(2->l) 


15m 




550 


4.7+0.3 


110+7 


1.6xl0 9 








HCN(1 


-0) 


30m 




550 


<0.33 


<2.0 


<l.lxl0 8 




PKS 0634-20 


220 


CO(1->0) 


12m 






<0.38 


<9.6 


<1.0xl0 9 


<1.6xl0 9 


B2 0648+27 


167 


CO(l- 


.0) 


12m 






<0.26 


<6.6 


<4.3xl0 8 


<6.4xl0 8 


B2 0722+30 


76 


CO(l- 


.0) 


12m 


0.0188 


450 


2.3+0.2 


58+4 


7.9x10 s 


1.2X10 9 






CO(2- 


►1) 


15m 




450 


6.8+0.2 


160+5 


5.4xl0 8 








HCN(1 




30m 




450 


<0.18 


<1.1 


<1.4xl0 7 




4C 29.30 


264 


CO(l- 


>0) 


12m 






<0.18 


<4.5 


<7.1xl0 8 


<l.lxl0 9 


Hydra A 


218 


CO(l- 


.0) 


12m 






<0.25 


<6.4 


<6.9xlO s 


<1.0xl0 9 


M 87 


19 


CO(l- 


.0) 


12m 






<0.46 


<11.7 


<1.0xl0 7 


<1.5xl0 7 


B2 1318+34 


94 


CO(l- 


.0) 


12m 


0.0233 


400 


7.0+0.3 


177+7 


3.6 xlO 9 


5.5xl0 9 


3C 285 


322 


CO(l- 


.0) 


12m 






<0.20 


<4.9 


<l.lxl0 9 


<1.7xl0 9 


PKS 1345+12 


503 


CO(l- 


.0) 


12m 


0.1220 


250 


0.24+0.06 


6.0+1.6 


3.3xl0 9 


4.9xl0 9 


PKS 1345+12 


503 


CO(l- 


>0) 


30m 




400 


1.6+0.08 


8.0+0.4 


4.3xl0 9 


6.5xl0 9 


3C 293 


182 


CO(l- 


.0) 


12m 


0.0448 


700 


2.2+0.09 


56+2 


4.3xl0 9 


6.4xl0 9 


OQ 208 


313 


CO(l- 


.0) 


12m 






<0.31 


<7.9 


<1.7xl0 9 


<2.6xl0 9 


3C 305 


166 


CO(l- 


.0) 


12m 






<0.34 


<8.6 


<5.4xlO s 


<8.2xl0 8 


TXS 1506+345 


181 


CO(l- 


.0) 


12m 


0.0449 


550 


2.5+0.09 


63+2 


4.8xl0 9 


7.1xl0 9 


PKS B1518+045 


208 


CO(l- 


.0) 


12m 






<0.45 


<11.3 


<l.lxl0 9 


<1.7xl0 9 


3C 321 


393 


CO(l- 


,0) 


12m 






<0.19 


<4.7 


<1.6xl0 9 


<2.4xl0 9 


3C 327 


429 


CO(l- 


.0) 


12m 






<0.22 


<5.4 


<2.2xl0 9 


<3.3xl0 9 


B2 1707+34 


329 


CO(l- 


.0) 


12m 






<0.19 


<4.9 


<1.2xl0 9 


<1.8xl0 9 


3C 390.3 


227 


CO(l- 


.0) 


12m 






<0.41 


<10.3 


<1.2xl0 9 


<1.8xl0 9 


3C 402 


102 


CO(l- 


.0) 


12m 






<0.41 


<10.3 


<2.5xl0 8 


<3.8xl0 8 


3C 403 


239 


CO(l- 


.0) 


12m 






<0.38 


<9.5 


<1.2xl0 9 


<1.9xl0 9 


3C 405 


225 


CO(l- 


>0) 


30m 






<0.37 


<1.9 


<2.1xlO s 


<3.2xl0 8 


3C 433 


416 


CO(l- 


.0) 


12m 






<0.16 


<3.9 


<1.5xl0 9 


<2.2xl0 9 



a Luminosity Distance, calculated assuming Ho = 75 km s 1 Mpc 1 , qo = 0.5, and A = 0. 

b For all radio galaxies with CO(1^0) non-detections, T m bAti = 3T m b(rms)AupwHM, where T m b(rms) is the root-mean-squared main- 
beam temperature and Aiip^HM = 500 km s _1 . 

c Calculatcd assuming a = 1.5Mq [K km s _1 pc 2 ] -1 . 



CO IN RADIO GALAXIES 



Table 4 

Estimation of M dust , M dyn , and Column Densities 





Tdust 


M dust 


M(H 2 )/M dust 


Rco 


M dyn 


M(H 2 )/M dyn 


Column Density 


Source 


(K) 


(M ) 




(pc) 


(M ) 




(cm- 2 ) 


3C 31 


33 


2.7 x 10 6 


160 


75 


4.2 x 10 9 


0.10 


1.6 x 10 24 


3C 84 


64 


1.7 x 10 6 


1084 


154 


2.2 x 10 9 


0.83 


1.6 x 10 24 


3C 120 


55 


2.0 x 10 6 


884 


112 


7.6 x 10 9 


0.23 


3.0 x 10 24 


PKS 0502-10 


46 


3.8 x 10 6 


1075 


185 


1.3 x 10 10 


0.32 


2.5 x 10 24 


B2 0722+30 


50 


2.3 x 10 6 


525 


105 


4.8 x 10 9 


0.24 


2.2 x 10 24 


B2 1318+34 


58 


9.4 x 10 6 


572 


222 


8.0 x 10 9 


0.68 


2.3 x 10 24 


PKS 1345+12 


91 


8.1 x 10 6 


802 


194 


7.0 x 10 9 


0.94 


3.6 x 10 24 


3C 293 


40 


3.6 x 10 6 


1807 


222 


2.5 x 10 10 


0.26 


2.8 x 10 24 


TXS 1506+345 


54 


8.4 x 10 6 


860 


227 


1.6 x 10 10 


0.46 


2.9 x 10 24 


Average 






860 ± 460 






0.45 ± 0.29 
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